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Abstract 

High by-pass turbofan engines have fewer fan blades and 
lower rotation speeds than their predecessors. Conse- 
quently, the noise suppression at the low frequency end 
of the noise spectra has become an increasing concern. 
This has led to a renewed emphasis on improving noise 
suppression efficiency of passive, duct liner treatments 
at the lower frequencies. For a variety of reasons, pas- 
sive liners are comprised of locally-reacting, resonant 
absorbers. One reason for this design choice is to sat- 
isfy operational and economic requirements. The sim- 
plest liner design consists of a single layer of honeycomb 
core sandwiched between a porous facesheet and an im- 
pervious backing plate. These resonant absorbing struc- 
tures are integrated into the nacelle wall and are very ef- 
ficient over a limited bandwidth centered on their reso- 
nance frequency. Increased noise suppression bandwidth 
and greater suppression at lower frequencies is typically 
achieved for conventional liners by increasing the liner 
depth and incorporating thin porous septa into the hon- 
eycomb core. However, constraints on liner depth in 
modern high by-pass engine nacelles severely limit the 
suppression bandwidth extension to lower frequencies. 
Also, current honeycomb core liners may not be suitable 
for irregular geometric volumes heretofore not consid- 
ered. 

It is of interest, therefore, to find ways to circumvent 
liner depth restrictions and resonator cavity shape con- 
straints. One way to increase effective liner depth is 
to skew the honeycomb core axis relative to the porous 
facesheet surface. Other possibilities are to alter res- 
onator cavity shape, e.g. high aspect ratio, narrow chan- 
nels that possibly include right angle bends, 180° chan- 
nel fold-backs, and splayed channel walls to conform to 
irregular geometric constraints. These possibilities con- 
stitute the practical motivation for expanding impedance 
modeling capability to include unconventional resonator 
orientations and shapes. 
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The work reported in this paper is in the nature of 
a progress report and is limited to examining the im- 
plications of resonator axis skew on the composite nor- 
mal incidence impedance of an array of resonator chan- 
nels. Specifically, experimental results are compared 
with a modified impedance prediction model for high- 
aspect-ratio, rectangular, resonator channels with vary- 
ing amounts of skew relative to the incident particle ve- 
locity. It is shown that for resonator channel widths of 
1 to 2 mm, aspect ratios of 25 to 50, and skew angles 
of zero to sixty degrees, the surface impedance of test 
models can be predicted with good accuracy. Predicted 
resistances and reactances are particularly well corre- 
lated through the first resonance and first anti-resonance 
for all six test models investigated. Beyond the first 
anti-resonance, the impedance prediction model loses the 
ability to predict details of resistance and reactance but 
still predicts the mean trends very well. 


Nomenclature 

c sound speed, m/s 

h depth of honeycomb core, mm 

L channel length, mm 

R a aspect ratio 

Symbols: 

4> channel skew angle, degrees 

p ambient density, kg/m 3 

9 normalized acoustic resistance, real com- 

ponent of C, 

% normalized acoustic reactance, imaginary 

component of C 

C 9 + z'%, normal incidence acoustic 

impedance, normalized by pc 
Abbreviations : 
cc circular channel 

LaRC Langley Research Center 

SLA steriolithography apparatus 

SPL sound pressure level, dB 

TMM Two-Microphone Method 

ZKTL Zwikker-Kosten Transmission Line code 

An e' m time convention is used throughout this paper. 
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Figure 1. Nacelle liner installation. 


Introduction 

Although passive, absorptive liner technology has 
reached a relatively mature state, it remains a critical 
component in the overall noise abatement methodology 
for commercial aircraft engines. Further improvements 
and innovative applications for passive liner technology 
are likely in the future because of new developments in 
materials science, fabrication processes, and computa- 
tional modeling. It is these enabling technologies that 
provide the impetus for the project described in this pa- 
per. 

The advent of high by-pass turbofan engines has stim- 
ulated the search for and development of innovative liner 
concepts. These engines generally produce more noise 
than their predecessors in the low frequency part of the 
noise spectrum that relates to community annoyance. 1 
The broadband contribution from wide chord fans makes 
this noise less dominated by discrete frequencies. To 
complicate matters further, the high by-pass engine na- 
celle can accomodate less liner length relative to the na- 
celle diameter, especially in the fan inlet, as depicted 
in figure 1. Also, the nacelle wall thickness (which 
translates into depth available for liners) does not in- 
crease proportionately with increasing fan diameter. This 
evolving engine-cycle trend requires increasing liner effi- 
ciency just to maintain the same noise reduction obtain- 
able with the lower by-pass engines of the last century. 
One way to cope with liner depth constraint is to distort, 
or otherwise modify, the conventional honeycomb core 
to lower the resonance frequency and thus enhance low 
frequency absorption efficiency. 

The semi-empirical character of current impedance 
prediction models may limit the impedance prediction 
accuracy of liner structures that depart from the conven- 
tional honeycomb-core architectures. The reason for this 
is that the conventional honeycomb-core liner is typically 
modeled as a one-dimensional, non-dissipative acoustic 
system. The honeycomb core contribution to the liner 
impedance is purely reactive. The facesheet contributes 
nearly all of the resistance, as well as a mass-like reac- 
tive component that increases linearly with frequency. 
This model is incapable of accounting for the liner be- 
havior in a relatively small frequency range centered on 
anti-resonance frequencies. Thus, for a complex liner 


Conventional liner Modified liner 

(honeycomb core) (skewed resonator channels) 



h(cosc|)) = 0.71h at 4> = 45° 

-> 30% reduction in depth for no 
reduction in resonance frequency 


Figure 2. Conventional and skewed resonator channels. 


geometry where dissipation may be present in the res- 
onator channels and a composite impedance is desired 
for a group of contiguous resonators with channel-to- 
channel depth variations, this simple model may lose ac- 
curacy. Yet, the accuracy of impedance prediction mod- 
els for such structures is crucial to the achievement of 
target impedance spectra. It is therefore of some impor- 
tance that the impedance contribution from unconven- 
tional resonator cavity architectures be accurately mod- 
eled. The simplest departure from conventional res- 
onator geometry is a skewing of the channel axis rela- 
tive to the porous surface normal. A simple example of 
this is illustrated in figure 2. In the sketch on the left, a 
conventional honeycomb core of depth h, with a thin re- 
sistive facesheet, is shown. In the sketch on the right, the 
resonator channels are skewed by the angle (|) = 45°, re- 
sulting in a depth reduction of about 30%, but no change 
in resonance frequency. In some applications, not indi- 
cated in the sketch, channel-to-channel depth variations 
and splayed channel walls may be necesary to exploit ir- 
regularly shaped, available space. Depending upon the 
wavelength, the resonances of splayed wall resonators 
may be shape dependent if non-Helmholtz resonances 
are of concern. It is imperative that these effects be un- 
derstood if the desired optimum impedance spectrum in- 
corporating such unconventional resonator geometry is 
to be achieved. 

This paper is a progress report on the first phase of 
an investigation that focuses on the effects of resonator 
axis skew on normal incidence impedance of high as- 
pect ratio rectangular channels. The skew is measured 
relative to the resonator axis when it is coincident with 
the incident sound particle velocity at the resonator face. 
To achieve a measurable resistance without installation 
of a thin resistive facesheet in the plane of the resonator 
face, an array of rectangular channels is used with chan- 
nel airway widths ranging from 1 mm to 2 mm. Thus, the 
resistive component of the impedance is supplied by cav- 
ity wall dissipation. Although resonator channel skew is 
the main focus of this phase of the investigation, it was 
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(a) Configuration Cl ; Straight section 


incident plane 




(b) Configuration C2; Straight section + straight 





(e) Configuration C5; 30 degree section 


Figure 3. Sketch of test models. 


necessary for the predictive model to also account for 
high-aspect-ratio rectangular channels and for channel- 
to-channel depth variation. The necessity of including 
channel-to-channel depth variation was driven by test 
model fabrication convenience as will be described forth- 
with. Effects of cavity wall splay will be investigated in 
a later phase of this project. 

Test Model Fabrication and Design 

The test models were fabricated by means of a steri- 
olithography apparatus (SLA) that works to an absolute 
precision of ±0.0635 mm. The apparatus provides for a 
numerically-positioned laser beam to map out, in succes- 
sive layers, the prescribed solid portions of the test model 
in a fluid resin bath. The precision exposure of the fluid 
resin to the laser beam causes a layer-by-layer fabrication 
of the test model. A great advantage of this technology 
is that it allows fabrication of intricate, unitary structures 
that require no assemblage of machined parts. 


Six test models were fabricated to study the effects of 
channel skew on the normal incidence impedance of ar- 
rays of contiguous resonators in the form of rigidly ter- 
minated rectangular channels. The resonator incident 
planes are co-planar with the 5 1 mm-square impedance 
measurement plane of the NASA Langley Research Cen- 
ter (LaRC) Vertical Impedance Tube. The design philos- 
ophy for the test models will be discussed with the aid 
of the side-view sketches of figure 3 and the photograph 
of figure 4. The sketch sequence in figure 3 also rep- 
resents the testing sequence. The straight channel test 
models illustrated in the sketches of figure 3-a,b are the 
baseline test cases which are used, in part, to validate 
the high aspect ratio modification of the impedance pre- 
diction model discussed in the "Impedance Prediction 
Model" section. The effect of channel skew was first 
studied with the aid of the two test models shown in fig- 
ure 3-c,d. These test models are designed to compare 
effects of skewing half the length (L = 85.7 mm) of the 
straight channel array (of length 2L) about the center- 
span point by 30 and 60 degrees to generate test models 
respectively. (Note that these test models have 17 and 
10 channels for 30 and 60 degrees of skew respectively.) 
This causes the straight and skewed sections to interface 
in a mitered joint, shortening and lengthening the chan- 
nels on either side of the pivot point in a symmetrical 
manner. Thus, the averaged length of the skewed chan- 
nels, taken together, does not change with skew. Also, 
the length of the central channel (or partition, as the case 
may be) is invariant. In this manner, the test models of 
figure 3-b,c,d are designed to isolate the skew effect from 
possible "contamination" by the classical "end effect" as- 
sociated with the area discontinuity at the incident plane. 
The test models (figure 3-a,e,f) can be thought of as a 
replication of the test models (figure 3-b,c,d) with the 
upper straight section (of length L) removed. These test 
models are designed to provide a validation of the predic- 
tion model when end effects and skew effects coalesce at 
the incident plane. 

From a fluid dynamic viewpoint, the effect of skew 
is to turn the normal incident particle velocity through 
the angle of skew, cf>. It is easy to visualize that when 
this occurs at the mitered joint of the straight and skewed 
sections of the test models in figure 3-c,d above, then the 
particle velocity in the skewed section must increase by 
the factor 1 / cos 4> because of the corresponding reduc- 
tion in channel width. Intuitively, one would expect a 
particle velocity increase to cause an increase in acous- 
tic resistance contributed by the skewed section. The 
test models that combine straight and skewed channels 
through the mitered interface are designed to provide the 
simplest posssible validation of the skew effect because 
the skew effect is separated geometrically from the end 
effect at the incident plane. 
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An important design requirement for the test models 
was to achieve a measurable resistance using the NASA 
LaRC Two-Microphone Method (TMM) impedance 
measurement methodology. This methodology achieves 
its maximum accuracy over a resistance range of about 
0.25pc to lOpc and over a frequency range of 0.5 kHz 
to 3.0 kHz. This measurement envelop, together with 
the fabrication precision limit indicated above, drove the 
baseline test models (straight channel sections in fig- 
ure 3-a,b to have 20 channels with the following pa- 
rameters: channel airway width, 1 mm; channel depth, 
85.7 mm; and channel partition thickness, 1.63 mm (the 
minimum channel partition thickness). 

The airway width of all the skew sections (figure 3- 
c,d,e,f) was maintained at 1 mm with a fabrication preci- 
sion of ±0.0635 mm. However, the number of channels 
in the 30 and 60 degree skew section (both when joined 
with a straight section and stand alone) was proportional 
to the baseline test model cross-section area projected 
through 30 and 60 degrees, respectively. This approx- 
imated to 17 and 10 channels in the 30 and 60 degree 
skewed sections, respectively. Once these geometric de- 
sign choices are made for the skew sections, the channel 
airway and partition widths are determined for the asso- 
ciated straight continuation. For the straight section con- 
tinuation from the 30 degree skew section, the straight 
channel airway width is 1. 14 mm and the partition width 
is 1.96 mm. For the straight section continuation from 
the 60 degree skew section, the airway width is 1.98 mm 
and the partition width is 3.45 mm. All test model chan- 
nel arrays were terminated by a 13 mm-thick end-cap 
with a thin film of 0.08 mm Mylar sandwiched between 
the end-cap and channel array termination to serve as 
a sealing gasket between channels. It is believed that 
this arrangement provided for very nearly zero acous- 
tic particle velocity at each channel termination and at 
least 20 dB of acoustic isolation between channels. All 
channels had the same span- wise dimension of 5 1 mm to 
give aspect ratios (ratio of channel span-wise dimension 
to channel airway width) ranging from about 25 to 50. 

Impedance Prediction Model 

The set of test models described above were designed 
to validate an impedance prediction model for small- 
channel width, high aspect ratio, rectangular channels 
that included channel axis skew with respect to the inci- 
dent wave vector. The test model design parameters were 
chosen to exploit the existing experience with the NASA 
LaRC Zwikker-Kosten Transmission Line (ZKTL) prop- 
agation code for small diameter tubes. This code is de- 
scribed in previous reports. 2,3 The code is based on a 
summary by Tijdeman 4 of the classical work initiated 
in the latter part of the 18th century concerning sound 
propagation in small diameter tubes. The ZKTL-based 



Figure 4. Photograph of straight & skewed channel test 
models. 

impedance prediction code 3 is empirically corrected for 
rectangular channels with aspect ratios of approximately 
9 by multiplying the normalized resistance component 
by the ratio ( R a + l)/2; i.e., 

_ R a + 1 

Sec — 2 '“'cc i L Acc 

where R a is the aspect ratio, and the subscript “cc” in- 
dicates quantities computed using the smallest cross- 
section dimension as an equivalent circular channel di- 
ameter. 

Physically, the Zwikker-Kosten analysis allows one to 
calculate, in closed form, the propagation constant and 
characteristic impedance of a small diameter channel in 
which the acoustic boundary layer dominates the parti- 
cle velocity profile. This causes significant dissipation 
of the acoustic wave to ensue. In practice, this happens 
for air-filled channels if the "effective" tube diameter is 
on the order of 1 mm or less. Thus, the present SLA 
fabricated test models meet that requirement wihout ex- 
cessively compromising channel width precision. 

At NASA LaRC, the ZKTL impedance prediction 
model has been validated for arrays of parallel tubes with 
"equivalent" diameters on the order of 1 mm. In the 
present case with aspect ratios of 25 to 50, a generalized 
radius is used based upon the analysis of Richards. 5 This 
analysis indicated that propagation through non-circular 
channels with high aspect ratio cross-sections can be ad- 
equately described by the Zwikker-Kosten analysis if a 
generalized tube radius is employed. This generalized 
radius is given by twice the channel cross-section area 
divided by the cross-section perimeter. For a rectangu- 
lar channel, the generalized radius is the channel cross 
section area divided by the cross section perimeter. This 
generalized radius modification is incorporated into the 
the NASA LaRC ZKTL impedance prediction model for 
high aspect ratio rectangular channels discussed in this 
paper. 
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Vertical Impedance Tube Apparatus 


Impedance measurements were conducted in the 
NASA LaRC Vertical Impedance Tube apparatus (VIT). 
The essential features of the apparatus are depicted in fig- 
ure 5. The vertical orientation is convenient for measur- 
ing the impedance of unconsolidated granular bulk ma- 
terials that would otherwise require restraining to avoid 
spilling out into the impedance tube. The impedance 
measurement methodology at NASA LaRC is a version 
of the switchable two-microphone-method (TMM) that 
employs a rotating plug to switch the microphones and 
thus "average out" differences between the data acquisi- 
tion channels. This circumvents a tedious calibraton pro- 
cedure, the details of which are explained elsewhere. 6 ' 7 
The methodology relies upon both the incident and re- 
flected waves to be devoid of any measurable contami- 
nanton from higher-order mode content at the measure- 
ment microphone locations (i.e. plane waves only). This 
requirement implies that any reflected wave nonunifor- 
mity generated at the test surface must have no measur- 
able influence on the total wave field at the measurement 
microphone locations. 

The impedance measurement plane for the VIT is lo- 
cated at its exit (see figure 5) and the test specimen sur- 
face is made coincident with this plane. A reference mi- 
crophone is located about 6.35 mm from this plane and 
provides a measure of the total sound field at the test 
specimen surface. Figure 6 shows a photograph of the 
VIT with one of the present test models (3-c) mounted. 
Note that the transparency of the apparatus allows the test 
model surface and environs to be viewed. This feature 
is regarded as valuable for including oscillatory flow- 
visualization should the need arise. 

Discussion of Results 

Figures 7 through 12 show comparisons of impedance 
spectra measurements and predicted values for the six 
test models. All measurements were conducted with 
bandwidth limited white noise over the frequency range 
of 0.5 to 3.5 kHz, in increments of 25 Hz. The overall 
SPL (measured over frequency range of 0.5 to 3.0 kHz) 
was set to 140 dB at the reference microphone loca- 
tion. The impedance model predictions are calculated 
at 100 Hz intervals. 

Figures 7 and 8 are the results for the straight, 20- 
channel, test models of length L and 2 L, respectively. 
The dimension L was set to 85.7 mm so that at least one 
anti-resonance is well defined for the shortest test model. 
That is clearly the case for configuration Cl in figure 7 
where the first resonance at about 1 .0 kHz and a second 
resonance at about 3.0 kHz bracket the anti -resonance at 
about 2.0 kHz. On the scale of this graph, resistance val- 
ues of a few tenths of one pc are not resolved well in 
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Figure 5. Vertically oriented impedance tube (VIT). 



Figure 6. Photograph of Vertical Impedance Tube. 
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the vicinity or resonance. Note, however, that predicted 
values of both resistance and reactance match the mea- 
sured values very well throughout the frequency range 
of 0.5 to 3.0 kHz. For the straight, 20-channel, 2 L length 
test model (configuration C2), the results are shown in 
figure 8 and the same remarks generally apply as for fig- 
ure 7. In this case, there are 3 resonances (0.5, 1.5 and 
2.5 kHz) and 3 anti-resonances in the frequency range 
of 0.5 to 3.0 kHz. Measured and predicted impedances 
are again in excellent agreement. These results validate 
the impedance prediction model for contiguous, high- 
aspect-ratio, rectangular channel test models. 

Figure 9 shows results when a skew angle of 30 de- 
grees is introduced into the 17-channel test model of 
length 2 L (sketch 3-c, configuration C3). This figure 
should be compared directly with figure 8. Generally, 
both predicted resistance and reactance are in excellent 
agreement with measured results. Note that there is a 
significant decrease (but well predicted) in the measured 
values of resistance at the third anti-resonance for the 30 
degree skew test model. The additional secondary peak 
in predicted resistance at about 3.3 kHz was not vali- 
dated by measurments in the 51 mm-square impedance 
tube, since the maximum measurement frequency was 
3.0 kHz. 

Figure 10 shows results when a skew angle of 60 de- 
grees is introduced into the 2 L length straight section 
(sketch 3-d, configuration C4). This figure should be 
compared directly with figure 8. Again, the measured 
and predicted impedance values are in very good agree- 
ment. There are significant changes generated by the 
introduction of 60 degrees of skew, and the prediction 
model does an excellent job of predicting the general na- 
ture of these changes. For instance, the first resonance 
moves from 500 Hz to about 600 Hz for both prediction 
and measurement. In addition, the first anti-resonance re- 
sistance values are reduced for both prediction and mea- 
surement. It is interesting to note that the third anti- 
resonance essentially disappears with the addition of 60 
degrees of skew, a prediction that is validated by the mea- 
surement. 

Figures 1 1 and 12 show predicted and measured values 
of impedance for the the two test models without lead- 
in straight sections, i.e. the incident sound is directly 
incident on the skewed channels. These figures should 
be compared directly with figure 7. Here again, pre- 
dicted and measured impedance values are in excellent 
agreement. Note the good agreement between the pre- 
dicted and measured resistances at the peak of the anti- 
resonance at about 20 pc for the 30 degrees of skew and 
about 27 pc for 60 degrees of skew. 


Concluding Remarks 

The following observations are drawn from the current 
investigation: 

1 . Generally, measured and predicted impedance val- 
ues are in excellent agreement for both straight and 
skewed-channel test models. 

2. The prediction model fails to resolve details of mea- 
sured resistance variability at the high frequencies. 

3. The results clearly show that resonator channel as- 
pect ratio and skew factors have an important ef- 
fect on on the nomal incidence impedance. Further- 
more, the influence of these parameters can be ade- 
quately modeled by simple modifications to a well- 
known model for propagation through small diam- 
eter, circular tubes. The aspect ratio factor is in- 
cluded by constructing a generalized radius equal 
to twice the channel cross section area divided by 
the channel perimeter. The angle of skew effect can 
be equated to a decrease in the channel open area 
ratio by the factor 1 / cost]). 
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Figure 7. Normalized impedance for Configuration Cl 
(straight 20-channel baseline - length, L). 



Frequency, kHz 


Figure 9. Normalized impedance for Configuration C3 
(straight 17-channel & 30-degree skew). 



Frequency, kHz 


Figure 8. Normalized impedance for Configuration C2 
(straight 20-channel, baseline - length, 2 L). 



Frequency, kHz 


Figure 10. Normalized impedance for Configuration C4 
(straight 17-channel & 60-degree skew). 
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Figure 1 1 . Normalized impedance for Configuration C5 
(isolated 17-channel at 30 degree skew). 



Figure 12. Normalized impedance for Configuration C6 
(isolated 10-channel at 60-degree skew). 
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